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We study the scalar and spinor perturbation, namely the Klein-Gordan and Dirac equations, in 
the Kerr-NUT space-time. The metric is invariant under the duality transformation involving the 
exchange of mass and NUT parameters on one hand and radial and angle coordinates on the other. 
We show that this invariance is also shared by the scalar and spinor perturbation equations. Further, 
by the duality transformation, one can go from the Kerr to the dual Kerr solution, and vice versa, 
and the same applies to the perturbation equations. In particular, it turns out that the potential 
barriers felt by the incoming scalar and spinor fields are higher for the dual Kerr than that for the 
Kerr. We also comment on existence of horizon and singularity. 
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> ■ I. INTRODUCTION 



o 



For asymptotically fiat axially symmetric stationary electro-vacuum space-time, it is well known that the Kerr 



family is unique [1]. On relaxing the condition of asymptotic flatness, there occurs the NUT generalization of it, the 
, Kerr-NUT family [2-4]. It has recently been shown [5] that this family is also unique for the space-time admitting 
O ■ separable (Hamilton- Jacobi and Klein-Gordon) equations of motion. When separability is implemented a priori, it 
^ determines certain form of the metric functions which are the functions of radial and angle coordinates. Then the 



Einstein-Maxwell equations for electro- vacuum space-time readily admit the most general solution which turns out to 
be the Kerr-NUT solution [4] establishing its uniqueness in a straightforward and direct manner. When the asymptotic 
flatness is imposed, the NUT parameter vanishes yielding the unique Kerr family. The Kerr family is thus included 
in the more general Kerr-NUT family which is the most general family for a metric admitting separable equations 
of motion. It would thus imply separability of the equations of motion for the Kerr family and it is indeed, as is 
well-known, so [1]. 

The most remarkable feature of this method of obtaining the most general solution is that it exposes the duality 
transformation which keeps the metric invariant. The duality here means exchange of the mass and NUT parameters 
on one hand and the radial and angular coordinates on the other [5,6]. Else the Kerr-NUT solution has been known for 
long [4], but this invariance has thus far remained unnoticed. Like mass is the measure of gravitational electric charge, 
NUT parameter is supposed to be the measure of magnetic charge [7]. The Kerr-NUT solution is therefore a truly 
gravitational dyon solution [6] and it represents gravitational field of a rotating body having both the gravitational 
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electric (mass) and magnetic (NUT) charges. By the duality transformation, one can obtain the space-time dual to the 
Kerr solution.^ There exists duality between the field of a rotating electric charge (M) and rotating magnetic charge 
(NUT parameter, I) [9] and vice-versa. Though there exists pure NUT solution [2] without the rotation parameter 

0, but for the duality transformation its presence is essential. That is, the duality between gravitational electric and 
magnetic charges could be exhibited only when the body is rotating. 

The Kerr family has been studied extensively for motion of test particles as well as for the scalar, spinor, vector 
and tensor perturbations of the metric. In this paper, we would like to study the scalar and spinor perturbations, 

1. e. the Klein-Gordon and Dirac equations in the Kerr-NUT geometry. As we said earlier that by a simple duality 
transformation one can switch over from the Kerr solution to the dual Kerr (where M = 0). Prom the properties 
of the Kerr metric, it is possible to derive the corresponding properties of the dual Kerr metric. For instance, the 
scalar and spinor perturbation equations of the Kerr geometry can be translated to the corresponding equations for 
the dual geometry by the duality transformation which takes the Kerr metric to the dual Kerr metric. However, the 
duality works only at the equation level but not at the solution level; viz the solution of the scalar equation in the 
Kerr case can not be taken over to the solution of the corresponding equation in the dual case by any such simple 
transformations. This happens because the duality involves mass and NUT parameters and their presence and absence 
in the equation changes the nature of solution. In the Kerr metric there is mass parameter but no NUT parameter. 
Under the duality transformation, mass goes to NUT parameter and so the dual metric has NUT parameter but no 
mass. Hence the solution which is obtained with mass parameter present will have different character and would not 
go over to the one obtained with mass zero and NUT non-zero. In the present case, this has become more non-trivial 
due to the particular selection of variable transformation which reduces the equations into the form of wave equations 
whose solution is known. 

The duality would have been very interesting, had it worked at the solution level as well. That unfortunately does 
not happen here trivially. Considering the Kerr solution as the field of rotating gravitational electric charge (M) and 
its dual as the field of rotating gravitational magnetic charge (/ - the NUT parameter) . Then it should be of interest 
to study the scalar and spinor perturbations of those two fields and compare their behaviour. This is exactly what 
we wish to do. We should however note that the dual solution is indeed an exact solution of the Einstein vacuum 
equation which is asymptotically non flat and its asymptotic flat limit is flat spacetime (i.e. if we impose the condition 
of asymptotic flatness, it reduces to flat spacetime). We shall probe this space-time with the scalar and spinor fields 
and compare their motion with the corresponding motion in the Kerr geometry. Our main aim here is to study the 
perturbations first in the general Kerr-NUT space-time and then specialize to the Kerr and the dual Kerr cases by 
employing the duality transformation. For the duality to work it is required that a ^ 0. Hence we cannot go to a = 
limit of the pure NUT solution. As we are not able to establish any duality relation between the solutions in the Kerr 
and dual Kerr cases, we will compare the solutions graphically. 

The paper is organized as follows. In §11, we shall recall the Kerr-NUT metric and its invariance under the duality 
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transformation. Next two sections III and IV would respectively be devoted to the Klein-Gordon and Dirac equations 
and their transformation under the duality. In §111, we will also present the numerical solution of scalar perturbation 
equation. Then in §V, we shall present and discuss the various numerical solutions of the Dirac spinors which would 
be followed by discussion of horizon and singularity in §VI, and we end with a discussion of the main results. 

II. THE KERR-NUT METRIC AND THE DUALITY 

Here we study the interaction of particles of spin zero and spin half with the Kcrr-NUT black hole. We are familiar 
with the Klein-Gordan and Dirac equations in flat space by which one can investigate the behaviour of spin zero and 
spin half particles. In the curved space-time the form of these equations are modified. The Dirac equation was written 
in the Kerr geometry by Tcukolsky [10] and was separated by Chandrasckhar [11] into the radial and angular parts 
(for a comprehensive discussion, sec [1]). Recently, there has been a fresh investigation of the radial Dirac equation 
in the Scliwarzscliild, R,eissner-Nordstrom and Kerr geometries [12 14]. 

In the linearized test field approximation, scalar, vector and tensor perturbations of the Kerr geometry has been 
studied by several authors (e.g. [1,10]). The master equations [10] governing these linear perturbations for integral spin 
(e.g., gravitational and electromagnetic) fields were solved numerically by Press & Teukolsky [15] and by Teukolsky & 
Press [16]. Particularly interesting is the fact that whereas gravitational and electromagnetic reflected wave from the 
black hole was found to be amplified for certain range of incoming frequencies, however Chandrasckhar [1] predicted 
that no such amplification should take place for the Dirac waves because of the very nature of the potential experienced 
by the incoming fields. 

The space-time metric for the Kerr-NUT geometry [5] is given by 

ds^ = {dt - Pdcj>f + ^ [{F + l^)dcl> - adtY + ^dr^ + p^dO^ (1) 

where, 

F = + a^, f/2 = - 2Mr + a'^ + Ql-P,P = asin^e - 2lcose, = + X^,X = I + acos0. (2) 

Here, I is the NUT parameter, a and are rotation and electric charge parameters of the black hole. The most 
remarkable property of this metric is that it is invariant under the duality transformation, M il,r iX [5,6]. That 
is, interchange of the gravitational electric (M) and magnetic [l) charge also requires the corresponding interchange 
of the radial and angle coordinates. For the former, radial is the responding coordinate, while for the latter, it is the 
angle coordinate. Further the transformation, M ^ zZ,r <-> zA will take the Kerr solution to the dual Kerr solution, 
and similarly vice-versa. 

Now following the Newman-Penrose formalism, we introduce null tetrads (/, n, fh, m) to satisfy orthogonality rela- 
tions, l.n = 1, m.m = — 1 and l.rh = n.fh = l.m = n.m = 0. Thus we write the basis vectors of null tetrad in terms 
of elements of the Kerr-NUT geometry as 

If, = ^(C/2, V,0,-C/2P), = ^(C/2,p^0,-C/2p), 

rriu, = -^-={iasin6,Q,—p^,—i{F-\-l^)sinO), fha = _ ^ ^ {—iasin6,0,—p^,i{F + l'^)sin6) (3) 
pV2 p*\/2 
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and 

F = ^{F + f, U\ 0, a), ni^ = ^{F + f, -^7^ 0, a), 

m'^ = ^-j={iPcosec9,Q,l,icosec6), fh^ = _ ^ ^ {—iPcosecO, 0,1,— icosecO), (4) 
pv2 p*v2 

where p = r + iX and jf = r — iX. 

In the next two sections, we study perturbations to the Kerr-NUT (including the Kerr as well as the dual Kerr) 
space-time with the Klein-Gordan and Dirac equations and the corresponding propagation of scalar and spinor waves 
respectively. So the perturbation can be expressed as a superposition of stationary waves with different modes given 
as e'('^*+"*'^) , where a is the frequency of the waves and m is the azimuthal quantum number. 



III. THE KLEIN-GORDAN EQUATION 



The general equation of scalar field (wave) in curved space-time can be written as 

1 



=5^ iV^g^-'d^^) -ml^ = 



(5) 



where nip is the mass of the scalar field. Now substituting = e»(o'*+™<^)$ in (5) we get 

+ ^ + ^iU^<, .) + - + + 2lU^cotecosece 

dr ' 



sin 

2 1 



C/2 



stn^ 



sin9 89 



{(^fisinO) + nipP^^ = 0. 



Further choosing $ = Ro{r)So{6) and separating the radial and angular parts of (6) we write 



d 



QQ2 + ""t^-QQ sin^Q 



AlamcotOcosecB + + m^{l + acosS 



i2 1 



sin^O 



5o = 0, 



(6) 



(7) 



d a2 



2 2 



+ ^(2Z^-g^ + 2Mr)-A? 



Ro = Q, 



(8) 



where Ai is the separation constant. 

It is interesting to note that, following [5], if we perform the duality transformation M ^ il, r ^ iX and Rq ^ Sq, 
and with suitable redefinition of a, Eqns. (7) and (8) remain invariant. Further, under the duality transformation, 
M ^ il, r ^ iX and Rq So, Eqns. (7) and (8) for the Kerr solution (with I = 0) will go over to that of the dual 
Kerr solution (with M = 0). Similarly, the reverse duality transformation, I —iM, r ^ iX and Rq ^ Sq, will bring 
the equations back to the Kerr geometry. 

Let us choose the transformation of independent variables r and 9 as 



1, fl-cosO 

y = 2^"' (tt^9 



^2 

z = — -{r++ r-)— + r+r-r. 



(9) 



where r± = M ± + P — — Ql and further 



Z = U'^Ro. 



(10) 
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Thus from (7), (8), (9) and (10) we get 



'0 



+ V^^o = 0, 



(11) 



dy2 




(12) 



where 



Vy = Xlsin^O - P'^a'^ + Alamcosd + mlsin'^6{l + acosOf' - m^, 



(13) 



4(r - M) 




Here Vy and Vz are the effective gravitational potentials for angular and radial motion respectively of the scalar field. 
At 6* = and tt and the corresponding y = -co and oo, Vy reduces to —{2la — m)^ and —{2la + rnf respectively. 
From (13) it is clear that Vy scales with I and a. Keeping the other parameters unchanged if I and/or a increases so 
does the potential. In the case of the radial potential, Vz, it falls off at large r but it diverges at the horizon [see Eqn. 
(14)]. Physically it means that the perturbation at large distance should not affect the space-time near the black hole. 
With the perturbation propagating towards the horizon, its effect becomes important. Here, note that in the natural 
units (G = c = ft, = 1), all the parameters are dimensionless. 

Equations (11) and (12) arc the simple one-dimensional wave equations which we solve numerically and depict in 
Figs. 1 and 2 for the Kerr-NUT, Kerr and dual Kerr space-times with a certain choice of parameters. In principle, 
one should first solve the angular equation with the eigenvalue, Ai. Then inserting that value of Ai into the radial 
equation, one should solve the radial equation. However, at the first instance we would like to get a qualitative 
feeling of the solution and hence we have set here Ai = 1 throughout. This is a deficiency of the solutions presented 
here which are not complete in this sense. However, they do provide some useful qualitative insight. In future we 
plan to evaluate Ai and get the complete solution. Here our main aim is just to give an indication of form of the 
possible solution. We use the Runge-Kutta method with the natural boundary conditions at infinity. We know that 
at infinity the potential barrier is either zero or constant. Therefore the wave form should be sinusoidal which defines 
the boundary condition at infinity. For the angular solution, first note that Vy is independent of the mass M and 
hence it is the same for the Kerr-NUT and the dual Kerr space-times. From the Figs. 1 and 2, the variance of 
solutions as a function of the rest mass of the incoming scalar as well as the space-time parameters, a, I and M, is 
very clear. Figure 1 shows that the amplitude of the solution is significantly high for the Kerr case. This is because 
the net gravitational effect is stronger for the Kerr case than the Kerr-NUT. It should be noted that the NUT and 
rotation parameters tend to oppose each other. Also the solution shows the existence of singularity in the space-time 
ai 6 ^ (]. On the other hand. Fig. 2 indicates an extension of the inner space-time region down to r = for the dual 
Kerr and thus the wavelength of scalar waves would therefore be wider than that for the Kerr and Kerr-NUT metrics. 
For the radial mode, larger is the wavelength for less massive perturbation. 
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FIG. 2. Radial Klein-Gordan solution with a = 0.4, for (a) Kerr-NUT, a = 0.998, I = 0.99, (b) Kerr, a = 0.998, (c) dual 
Kerr, a = I = 0.998. Solid and dotted curves indicate the cases, nip = 0.4, 0.1 respectively. M = 1 for (a) and (b). 

Note that the Klein- Gordon equation remains invariant under the duality transformation and the equation for the 
Kerr geometry can be transformed to that for the dual Kerr geometry and vice versa, but as is mentioned in §1 there 
is no easy way to exploit explicitly this duality for the solutions. One of the reason behind this is due to the certain 
choice of the independent variables to cast the equation in the familiar simple form which facilitates the analysis of 
motion through the behaviour of potential. The angular potential, Vy, depends upon the NUT parameter, I, and is 
free of the mass of the black hole. The Kerr solution has only kinematic part coming from the rotation parameter, a. 
While the radial potential, Vz, depends upon both M and I. We can go to the corresponding potential and solutions 
for the Kerr and the dual Kerr space-times by setting I = and M = respectively in the Eqns. (13) and (14). 

IV. THE DIRAC EQUATION 



Following [14] , derivative of spinor can be written as 

D^P^ = d^P^ + iqA^P^ + T^,P\ 



(15) 



where and T'f^^ are electromagnetic and gravitational gauge fields respectively. Thus, the Dirac equation in the 



Newman-Penrose formalism can be written as 



a'Xs'D^P^ + iHpQ^ ec'B' = 0, (16) 



^ab'-OmQ"" + il^pP^'ecB' = 0. (17) 

For a vector Xi, a\g,Xi = Xab'',A, B = 0,1; and is the mass of the Dirac particle. In terms of this new basis 
of the Newman-Penrose formalism, the Pauli matrices can be written as 



(18) 



Following [1] and writing 

P° = i^i,pi=F2,Q'' =Gi,Q° = -G2 

we get 

F{d^ + iqA^)Fi + m''(5^ + iqA^)F2 + (e - p)Fi + (tt - a)F2 = inpGi, (19) 
m'^id^ + iqA^)Fi + 71^(9^ + iqA^)F2 + {fi - 7)i=2 + {P - t)Fi = i/XpGa, (20) 
Fid^ + iqA^)G2 - m'*(5^ + iqA^)Gi + {e* - p*)G2 - (tt* - a*)Gi = i/XpFa, (21) 

n^{d^ + iqA^)Gi - in^'{d^ + iqA^)G2 + (m* - 7*)Gi - (/3* - t*)G2 = iHpF^. (22) 

These arc the Dirac equations in the curved space-time. 
We define the spinor fields as follows: 

and the electromagnetic vector potential, = {At,0,0,A^). For the Kerr-NUT metric, it can be expressed as 

A, = A, = (24) 

We also list the spin coefficients in terms of the Kerr-NUT metric as 

1 cot9 iasinO 
5= , 3 = — TT = —p= , e = 0, 

iasinO r — M _ . ^. 

Now by substituting Eqns. (23-25) into Eqns. (19-22), we reduce the Dirac equation to 

Do/i + 2-5£i /2 = [ippr + pp{l + acose)]gi (26) 
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where, 



and 



U^V\ f2 - 25£t = -2[ifipr + Hp{l + acose)]g2 (27) 

2 2 

Vog2 - 2-ic\gi = [i/ipr - fj,p{l + acose)]f2 (28) 

2 

U^V\ gi + 25£ 1 = -2[i/ipr - iJ.p(l + acose)]fi (29) 

2 2 

I^J, = |: - ^ (^(i^ + + am + gg.r) + 2n^l^ (30) 



= — + (crP + m)cosecO H- ncotO, 

fin 



de 

£t = A _ ((jP + m)cosece + ncotO. (31) 

Further writing /i(r,0) = R_.{r)S_._{e), h{r,0) = Riir)Si{0), gi{r,e) = Ri{r)S_i{0), g2{r,e) = 
-R_i (r)5_|_i (0), we can separate Eqns (26)-(29) into angular and radial parts as 

'Ci/2'S'i/2 = -(A2 - mp{l + acos6))S_i/2, £1^2*^-1/2 = (-^2 + mp{l + acos9))Si/2, (32) 
U'DoR_i/2 = {^2 + irUpT 

)URi/2, UVlURy2 = (A2 — impr)R_i/2- (33) 

Here, Trip is the normalized rest mass of the incoming particle and A2 is the separation constant. Unlike the case of 
the scalar perturbation, the invariance of the spinor perturbation equations under the duality transformation is non 
trivial, as the up and down spinor components couple with each other. However, if we choose the duality relations 
between spinor components as iS'1/2 ^ ^^ ""^^^-^-1/2 ^-nd 5'_i/2 ^ U^^^Ri/2, Eqns. (32-33) do remain invariant when 
M ^ il, r iX. As before, then we can switch over the Dirac equation in the Kerr space-time to that in the dual 
Kerr space-time and vice versa by the duality transformation. Further the interesting point to note is that under the 
duality transformation, the up-spinor transforms to the down and vice versa, with the multiplication/division of U. 
Next we study the separated angular and radial equations. 

IV. A. Angular equations 



If we decouple the equations in (32) for ^.i, we get 

■m„asin9 



(34) 



r 1 rt 5_ 1 + "V S_ 1 + (Ai - ml{l + acosef)S_ 1 = 0. 

2 2 2 \ + mp{l + acosO) 22^^ ' ' 2 

Similarly one can decouple the equations in (32) for 5*1 . Let us choose 

= (X2 ± rupl ± ^) log{l - cosO) - (\2 ± nipl T ^) log{l + cosO). (35) 
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Thus the decoupled equations for S_i [i.e. Eqn. (34)], and for Si [not given here], can be reduced to 

+ 2 



-^+W^S,.=0 (36) 
where 

cot'^9 



W-. 



±{aP + m)cot9cosec9 + acosec6{2asin6cos0 + 2lsin6) h 



^ {\2 ± mp{l + acosO)) 

/ ^ ^o .0 o,, m„asin0 ( cotQ , „ , ^\ 

- (crP + mYcosec^e + Xi- miU. + 00056*)^ ± - — ■ -r — — T (o-P + m)cosece 

A2 ± mp(< + acosC) \ 2 / 



2 4 

(37) 



Eqn. (36) is a simple one dimensional wave equation. Here again the angular potential depends upon / and is free 
off M. 

The above gravitational potentials are related to the physical parameters and variable in a rather complicated 
manner compared to the case of scalar perturbation. For both ^ — > and tt, angular potential barriers reduce to zero 
(W± — > 0). We can study the potential barriers felt by the incoming spinors with different masses and frequencies. In 
principle one can choose any mass and frequency of the field but to bring a significant interaction with the black hole 
we choose the mass and frequency in such a manner that the Compton wavelength of the incoming fields becomes 
same order as the radius of the outer horizon to the black hole. Also the frequency of the perturbation should be 
same order as the inverse of the light crossing time to the radius of the black hole. Thus we choose 

a ~ rup ~ r^^. (38) 

Figures 3 and 4 indicate the behaviour of potentials felt by the spinors in the angular direction in the transformed 
coordinate system. For a particular a, as nip decreases the coupling strength between the spinor and space-time 
curvature reduces then the peak of the barrier goes down. Similarly for the same rUp if a reduces the peak goes down. 
The interesting feature that emerges out is that the height of the barrier is quite sensitive to the spin orientation 
(up/down) of the incoming spinor field. This is in accordance with the expectation that a spinning black hole (with 
a particular spin orientation) can distinguish the up and down spin of the perturbation separately. Further a spinor 
feels higher potential barrier when it is aligned with the black hole. We also note the mirror inversion symmetry 
between W+ and W-. 
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FIG. 3. Variation of W+ as a function of m+, M = 1 when (a) a = 0.998, I = 0.99, (b) a = 0.998, / = 0.1, (c) a = I = 0.5, 
(d) a = 0.5, / = 0.1. Solid curves indicate the cases, a = 0.4; nip = 0.4,0.3,0.2,0.1,0 from top to bottom and dotted curves 
indicate a = 0.1; rup = 0.1, from top to bottom. 
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FIG. 4. Variation of W- as a function of M = 1 when (a) a = 0.998, I = 0.99, (b) a = 0.998, I = 0.1, (c) a ^ I = 0.5, 
(d) a = 0.5, / = 0.1. Solid curves indicate the cases, a = 0.4; nip = 0.4,0.3,0.2,0.1,0 from top to bottom and dotted curves 
indicate a = 0.1; rup = 0.1, from top to bottom. 



IV. B. Radial equations 



Decoupling the equations in Eqn. (33) for R_i we write 



U^V\ V„R_ 1 - JlIh!E^VoR_ 1 - (A^ + my)R_ i = 0. (39) 



However, following the earlier works [1,13], we cast this into the form. 



which is easier to attack. Here we define 



+ (ct^ - V±) Z±=0 (40) 



Z± = C/i?le*e/2^i^_le-^®/^ e=tan-i(^), (41) 

2Mr+-Q2+2Z2 + 2I2±92iI+ fr \ 2Mr- -Ql + 21^ + / 

r* = r -I log 1 log 

r+ — r_ \r+ J r+ — r- \ r_ 

1 1 /"mj,r\ 
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[a;2(Ai + m2r2) + ^i^]2^ V2 p ; A2 p ; p si 



^ [a;2(Ai + m2r2) + ^-^p 



2r H I (A2 + m^r ) + 2mpU) r H ^-^ 



a 



(43) 



where 



a;2=r2 + a^ = + f + ^ + SR^. (44) 

(T (T 

From Eqn. (42), it is clear that for 

qQ*r++am 

coefficient of log ^ vanishes and thus for a < ag, the relation ft: — r becomes multivalued. Prom Eqn. (43) it 
is clear that the potential, V±, diverges at r = \a\. This regime where V± diverges, is well known as super-radiance. 
For a > as, \a\ < r+, and the wave can not reach at r = |a| and hence there is no question of super-radiance. For 
a = as, \a\ = r+, and the potential [Eqn. (43)] diverges exactly at the event horizon. For a > as, f* runs from —00 
to 00 as r varies from r+ to 00; but for a < as, ad, both r ^ r+ and r ^ 00, r* ^ 00. The turning point of the f^—r 
relation occurs at r = \a\. Thus the analogue of ergo-sphere is defined by the region r+ < r < \a\. However, the 
spinor particles do not experience super-radiance and hence no energy can be extracted by the spinor perturbation 
(e.g. [1,13,18]). However the result is non null for the electromagnetic perturbation and energy is possible to extract 
from the black hole. 

It is clear from Eqn. (43) that, like angular potential barrier, here also the gravitational potentials are related to the 
physical parameters and variable in a complicated manner. Close to the horizon, the potential barrier, V±, reduces 
to zero but at a large distance V± — > m^. Thus from Eqn. (40), it appears that a has to be greater than or equal to 
nip, otherwise spinors can not enter into the gravitational field. 

Figures 5 and 6 show the behaviour of potential barriers (V±) for diff'erent choices of the frequency and mass of 
the spinor. Here also we choose a rUp ~ r^^, to consider a significant interaction between the incoming spinor and 
the black hole. For a particular frequency, as the mass of the spinor decreases the asymptotic height of the barrier 
decreases. 
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FIG. 6. Variation of V- as a function of f, for a > as, when (a) a = 0.998, I = 0.99, (b) a — 0.998, / = 0.1. Solid curves are 
for a = 0.7 and dotted curves for a = 0.4. From the top to bottom rup runs through 0.4, 0.3, 0.2, 0.1, 0, and M = 1, m = —1/2, 
Q, =0. 

In Fig. 7, the behaviour of potentials (V+) in super-radiance regime is shown. Clearly, at a certain r the potential 
diverges. For a particular a, if rUp increases, the singular point shifts outward. The behaviour of corresponding V- is 
similar. 
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FIG. 7. Variation of V+ as a function of f, for cr < cts (cr = 0.1), when (a) a = 0.998, / = 0.5, (b) a = 0.998, I = 0.1. Solid 
and dotted curves indicate lUp = 0,0.1 respectively with M = 1, m = —1/2, Q» = 



IV. C. Comparison of Kerr with dual Kerr 
For the dual Kerr solution M = 0, the corresponding f* would read as 




for a = l, and 

for a > Z. The radial potential would depend upon a and I only, keeping unchanged the form of Eqn. (43). 

Figures 8 and 9 show the variation of angular potential in the Kerr and dual of Kerr metrics. It is seen that the 
peak of the barrier (and the maximum height of it) in the dual Kerr for a particular set of a and nip is higher than 
that in the Kerr case. However W± in the dual Kerr is the same as for the Kerr-NUT for the angular motion because 
of its M independence. As before, here again W+ is higher than W- because of the opposite spin orientation of the 
perturbation, particularly in the higher frequency regime. 
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FIG. 8. Variation of W+ for the dual Kerr (M = 0) case as a function of u+ for (a) a — 0.998, I = 0, M = 1 (Kerr) 
(b) a = I = 0.998. From top to bottom, solid curves are for a = 0.4; nip = 0.4,0.3,0.2,0.1,0 while dotted are for a = 0.1; 
rup = 0.1, 0, and m = — 1/2, Q* =0. 
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FIG. 9. Potential is plotted against U- for (a) o = 0.998, / = 0, M = 1 (Kerr) (b) a = Z = 0.998, M = (dual 
Kerr). From top to bottom, solid curves are for o = 0.4; rup = 0.4, 0.3, 0.2, 0.1, and dotted are for a = 0.1; rup = 0.1, 0, and 
m = -1/2, Q* =0. 



In Figs. 10 and 11, we compare the radial gravitational potential in the Kerr and dual Kerr cases. It is again seen 
that the peak of the barrier (and the maximum height of it) in the dual Kerr for a particular set of a and rUp is 
higher than that in the Kerr. In the dual Kerr, for a particular rUp, if a decreases, the peak of the barrier distinctly 
increases. 
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FIG. 11. Variation of V- as a function of f., when (a) a = 0.998, / = 0, M = 1 (Kerr) (b) a ^ I = 0.998, M = (dual 
Kerr). Solid curves are from top to bottom for a = 0.7; rrip = 0.7,0.6,0.5,0.4,0.3,0.2,0.1,0 and dotted are for a = 0.4; 
rup = 0.4, 0.3, 0.2, 0.1, 0, and m = -1/2, Q« = 0. 

Therefore, from the above discussion we now have a clear picture of the behaviour of gravitational potential barriers 
felt by the spinors in the Kerr-NUT, Kerr and dual Kerr space-times. We will next study their solutions. 

V. SOLUTION FOR SPINOR PERTURBATION 

Here we like to solve Eqns. (36) and (40) for the various sets of the physical parameters. In view of the stationary 
and axisymmctric nature of the background space-time, it is natural to write the perturbation as superposition of 
waves of different modes, exp[i{(Tt + rn^i)]. Wc have thus to solve the radial and angular equations. These equations 
are however not solvable analytically and hence we arc forced to seek for numerical solutions. We shall employ the 
well known Rungc-Kutta method with the boundary conditions at infinity. With the aim of getting the qualitative 
feeling of the solution, wc shall again set A2 = 1, which should in principle be evaluated in an exact manner. The 
complete sohition has been obtained for the Kerr geometry [17]. We shall however defer this for a future consideration 
for the Kerr-NUT geometry. 
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The numerical solutions are shown in the following figures. The boundary condition at infinity is fixed by demanding 
sinusoidal wave with the wave number, k± = (j^ — V±. 

Figure 12 represents the behaviour of spinors in the angular direction for the Kerr-NUT and Kerr geometry. Recall 
from Figs. 3 and 4 that potential (which is also the square of wave number according to our notation) is positive only 
in the vicinity of ^ = 7r/2, then (36) can admit a harmonic oscillator like solution. We will essentially be interested in 
this region only because the regions of negative W± can only give exponentially increasing and/or decreasing solutions. 
The solutions clearly bring out the role played by the NUT parameter. 
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FIG. 12. Angular spin-up and spin-down solutions with a = nip = 0.4, for (a) Kerr-NUT, a = 0.998, I = 0.99 (b) Kerr, 
a = 0.998. (c) Kerr-NUT, a = I = 0.5, (d) Kerr, a = 0.5. Solid and dotted curves indicate the up and down spinors respectively 
as marked in the figure too. Other parameters are m = — 1/2, Q* =0. 



Figures 13 and 14 depict the corresponding radial evolution of spinors. The opposite sign of spin-spin coupling 
between the black hole and spinors of up and down alignment reflects into the solutions which show the amplitude of 
the spin down wave is significantly less compared to that of the spin up. The solutions in the dual Kerr (Fig. 13c) 
are shown in the semi-log scale just to catch the abrupt variations of curves as r — > in the same plot extended upto 
a large r. At r — > oo, V± — > m^, therefore the wave number of spinors becomes small at far away from the event 
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horizon which produces a large wavelength compared to that at a small r. For a = 0.4, cases with a = 0.5 lie in the 
super-radiance regime for the dual Kerr and then the corresponding radial potentials diverge. 
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FIG. 13. Radial spin-up and spin-down solutions with a = nip — 0.4, for (a) Kerr-NUT, a = 0.998, I = 0.99 (b) Kerr, 
a = 0.998, (c) dual Kerr, a = I — 0.998. Solid and dotted curves indicate the up and down spinors respectively as marked in 
the figure too. Other parameters are m = —1/2, Q* = 0; for (a) and (b) M = 1. 
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FIG. 14. Radial spin-up and spin-down solutions with a = rUp = 0.4, for (a) Kcrr-NUT, a — I = 0.5, (b) Kerr, a = 0.5. Solid 
and dotted curves indicate the up and down spinors respectively as marked in the figure too. Other parameters are m = —1/2, 
Q, =0, M = l. 

Now we will compare the solutions of the spinor perturbation equations in the various regimes of incoming spinors 
so that the black hole will appear to act as a mass-spectrograph. Figures 15 and 16 bring this feature for the angular 
up and down spinors respectively. As before only that range of 6 is chosen when the solution is harmonic and it 
however covers most of the range of 9. The deviation of solutions for different sets of {a, nip} is very clearly shown in 
the figures for the Kerr and Kerr-NUT space-times. 
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FIG. 16. Angular spin-down solutions with o = 0.4, rup = 0.1 (solid curve); cr = 0.4, = 0.4 (dotted curve); a — 0.7, 
rup = 0.4 (dashed curve); for (a) Kerr-NUT, a = 0.998, / = 0.99, (b) Kerr, a = 0.998, (c) Kerr-NUT, a = l = 0.5, (d) Kerr, 
a = 0.5. Other parameters are m = —1/2, Q* = 0. 

Figures 17 and 18 indicate the solutions for the up and down spinors respectively in the radial direction in case of 
the extreme Kerr-NUT, Kerr and dual Kerr space-times for different choices of {a,mp}. The most interesting fact 
in all the figures is that when a ~ rUp (dotted curves) the wavelength as well as the amplitude of the spinors are 
high compared to two other cases. This is physically understood as a rUp indicates the very small wave number 
at least at a large distance, that trend continues up to the inner radii, except at very close to the horizon where the 
strong effect of a black hole suppresses any individual identity of the matter. For a > rUp, the perturbation of higher 
frequency gives rise to the solution of lower wavelength which is physically understood. Again the solutions for the 
dual Kerr are distinctly different from others. As before the overall amplitudes of up spinors which are oriented with 
the black hole's spin in a parallel sense, are high compared to that of down spinors oriented in an anti-parallel sense. 
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FIG. 17. Radial spin-up solutions with a — 0.4, nip — 0.1 (solid curve); a = 0.4, nip = 0.4 (dotted curve); a = 0.7, mp = 0.4 
(dashed curve); for (a) Kerr-NUT, a = 0.998, / = 0.99 (b) Kerr, a = 0.998, (c) dual Kerr, a = l = 0.998. Other parameters are 
m = -1/2, Q* = 0; for (a) and (b) M = 1. 
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FIG. 18. Radial spin-down solutions with a = 0.4, = 0.1 (solid curve); a = 0.4, rUp = 0.4 (dotted curve); a = 0.7, = 0.4 
(dashed curve); for (a) Kerr-NUT, o = 0.998, / = 0.99 (b) Kerr, o = 0.998, (c) dual Kerr, a = l = 0.998. Other parameters are 
m = -1/2, Q* = 0; for (a) and (b) M = 1. 

The similar sets of solution are compared in Figs. 19 and 20 but for an intermediate rotating black hole. All 
the basic features are similar to that of Figs. 17 and 18 which are not repeated again. For the black hole angular 
momentum (Kerr parameter) a = 0.5 the corresponding as = 0.5 in case of the dual Kerr with I = 0.5. Therefore 
only the a = 0.7 case is out of the super-radiance regime and depicted in Figs. 19c and 20c. 
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FIG. 19. Radial spin-up solutions with a — 0.4, trip = 0.1 (solid curve); a = 0.4, nip = 0.4 (dotted curve); a = 0.7, rUp = 0.4 
(dashed curve); for (a) Kerr-NUT, a = I = 0.5 (b) Kerr, a = 0.5, (c) dual Kerr, a = I = 0.5. Other parameters are m = —1/2, 
= 0; for (a) and (b) M = 1. 
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FIG. 20. Radial spin-down solutions with a = 0.4, nip = 0.1 (solid curve); a = 0.4, nip — 0.4 (dotted curve); a = 0.7, nip = 0.4 
(dashed curve); for (a) Kerr-NUT, a = I = 0.5 (b) Kerr, a = 0.5, (c) dual Kerr, a = I = 0.5. Other parameters are m = —1/2, 
Q, = 0; for (a) and (b) M = 1. 



VI. HORIZON AND SINGULARITY 



The distinguishing feature of the Kerr-NUT geometry is the presence of the magnetic charge which is called the NUT 
parameter, /, that causes the space-time to be asymptotically non-flat. The horizons are defined by C/ = and hence 
they will occur at r± = M± ^JM"^ — o? — Q^. We shall now switch off the electric charge, i.e. Q* = 0, as it is not 
very pertinent to carry it through for the purpose in question. The singularity would occur where curvatures diverge 
and that would happen when acosfff' = 0- That means singularity would be located at r = 0, cos6 = —I /a. 

Clearly there occurs no singularity when I > a. It is the ring like (r = 0, ^ = 7r/2) for the Kerr case when I = and 
string like when a > I > (Kerr-NUT case). The string like character of the singularity is the characteristic of the 
presence of the NUT parameter [19]. 

On the other hand, horizon would always occur for + P > a^. For I > a, there occurs an interesting situation 
of having a horizon without singularity. The space-time is then non singular and regular everywhere. This is very 
strange because horizon should generally cover singularity. It has happened due to dominance of the NUT parameter 
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over the Kerr parameter. It is however well known that the NUT parameter also gives rise to the closed time-like 
geodesies. This suggests that the NUT parameter should never dominate over the Kerr parameter to ensure physical 
reasonableness. 

If however, we wish to have the canonical picture of a black hole space-time having the horizon covering a singularity, 
then P < a^, which would imply a ^ always. Recall that this was also required for the duality transformation 
between the electric (mass) and magnetic (NUT) charge. In case of the dual Kerr solution M = and the horizon 
would be given by r = i-v/Z^ — which would require > a^, while existence of singularity would require l"^ < c? . 
Requiring both the horizon and singularity to coexist, we are unambiguously led to = c? . The horizon occurs at 
r = and it becomes singular when ^ = 0, tt. The dual Kerr solution may not be very realistic but it provides an 
interesting vacuum space-time which is free of the usual electric gravitational charge (mass). The most remarkable 
feature is that general relativity admits a truly gravitational dyon solution which has the magnetic limit as the dual 
Kerr solution. 

VII. DISCUSSION 

We have studied the scalar and spinor perturbation in the Kcrr-NUT (which includes the Kerr and dual Kerr) 
space-time which describe gravitational field of a rotating gravitational dyon with electric and magnetic charge. It is 
known that the Kcrr-NUT space-time is invariant under the duality transformation which exchanges mass and NUT 
parameters as well as radial and angle coordinates [5,6]. By the duality transformation one can go from the Kerr to 
the dual Kerr solution and the vice versa. Interestingly, this duality transformation requires the rotation parameter 
to be non zero. Here we have shown that the same is true for the Klein-Gordon and Dirac equations in the Kerr-NUT 
space-time. That is, they are invariant under the duality, and they transform from the Kerr background to the dual 
Kerr, and the vice versa under the duality transformation with some appropriate rescaling of parameters. 

In the Kerr geometry, angular part of the equation is free from mass and simply involves the kinematic aspects 
arising from rotation. With the NUT parameter, angular part attains active dynamical meaning by its presence. 
The radial part would however involve both the mass and NUT parameters. The behaviour of potentials and their 
sohitions in different cases are shown in various figures. Potential barriers are higher for the dual Kerr space-time. 
This is because of the absence of mass, which produces the usual l/r attractive potential while the NUT parameter, 
Z, contributes asymptotically. Thus gravity is as expected stronger for the Kerr than the dual Kerr. 

Though the Klein-Gordon and Dirac equations could be transformed under the duality transformation from the Kerr 
to the dual Kerr case, their solutions could not be transformed into each other by any simple duality transformation. 
Had the duality worked for the solutions of the equations, it would have indicated something profound. The problem 
is that in the Kerr background I = while for the dual Kerr M = 0, the character of the equation changes drastically 
in the two cases and hence their solutions can not be related by the duality transformation. This however does not 
rule out the possibility that there may occur some other transformation which may work with the solutions. Note 
that in obtaining the solutions in a known form, we have defined the new independent variables in a complicated form 
which does not let duality transformation work at the solution level. With a proper and suitable new definition, it may 
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be possible to relate the solutions by the duality transformation. This as well as finding the complete solution with 
Ai , A2 properly determined (rather than the qualitative solutions as considered here) of the perturbation equation we 
leave for future study. But we do however believe that this would not significantly alter the qualitative character of 
the solutions. 
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